Potential energy landscape-based extended van der Waals equation 
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The inherent structures (IS) are the local minima of the potential energy surface or landscape, 
U(r), of an N atom system. Stillinger has given an exact IS formulation of thermodynamics. Here 
the implications for the equation of state are investigated. It is shown that the van der Waals (vdW) 
equation, with density-dependent a and b coefficients, holds on the high-temperature plateau of the 
averaged IS energy. However, an additional "landscape" contribution to the pressure is found at 
lower T. The resulting extended udWequation, unlike the original, is capable of yielding a water-like 
density anomaly, flat isotherms in the coexistence region vs vdW loops, and several other desirable 
features. The plateau energy, the width of the distribution of IS, and the "top of the landscape" 
temperature are simulated over a broad reduced density range, 2.0 > p > 0.20, in the Lennard- Jones 
fluid. Fits to the data yield an explicit equation of state, which is argued to be useful at high density; 
it nevertheless reproduces the known values of a and b at the critical point. 



I. INTRODUCTION 

Classical thermodynamics and dynamics are ulti- 
mately governed by the potential energy surface or land- 
scape, denoted PEL, denned at constant volume in the 
3iV dimensional configuration space of all the atomic co- 
ordinates of an atomic or molecular system. The canon- 
ical partition function is determined by the integral of 
the Boltzmann factor over the configuration space. Still- 
inger and Weber p| proposed a novel approach to this 
longstanding theoretical challenge. The space is parti- 
tioned into the basins of attraction of the local minima, 
named 0] inherent structures (IS), and the integral be- 
comes a sum over basins. Knowledge of distribution in 
energy of distinguishable IS and of the "vibrational free 
energy of the basins, A v (U is ,T, V) then allows 0,0,13 
evaluation of the partition function, the Helmholtz free 
energy A(T, V), and all of thermodynamics. Of course, 
obtaining these PEL properties is difficult, and detailed 
implementation of the IS formalism is in its infancy. 

Intensive quantities are defined 0, 0] as 4>=Ui S /N , 
a v =A v /N. The IS distribution is £l(</>) = Cexp(Na(4>)), 
where C has dimensions of inverse energy. Then, 

A(Ndp, T, V)/N = 4> + a v (<t>, T, V) - k B Ta(cp, V), (1) 

and a plausible, but not unique 0, expression for 
the configuration entropy/particle is determined by a, 
S c (N(j),V)/N = k B cr(<f>,V) = s c {<j>,V). In the limit of 
large TV the IS energy which minimizes A, the thermody- 
namic average denoted <j)(T,V), is overwhelmingly dom- 
inant. The thermodynamic A is obtained by replacing <p 
with <f>(T, V) in EqEl 

The equation of state 0, is determined by the rela- 
tion, P = -(dA/dV) T = +p 2 (d(A/N)/dp) T , where p is 
the number density; we will use density instead of volume 
in the following. Given <p(T,p), a((j),p), and a v (<j),T, p), 
one can simply evaluate Eq ?? and differentiate. Alter- 
natively, note that the density derivative of A acts both 
on the explicit p-dependence of a and a v , and on that 
implicit in 4>(T,p). Using (dA/ d<p)=Q at equilibrium it 



may be seen that the implicit contributions sum to zero 
and, 



p = p 2 {{^mT,P),T, P ) 

dp 

+(g)(#T,p),I», 



(2) 



including the explicit p-dependence only. 

Little is currently known about these quantities. A 
Gaussian approximation, 



a(^,p) = a-(^-Mp)) 2 /^ 2 (p), 



(3) 



where <^o is the band center and S 2 is the standard devia- 
tion squared/atom, is reasonable for fluids if one believes 
|E 0> IS 13 that the IS are built up from weakly interacting 
local regions. The total number of IS is p] oc exp(aN) . 

A harmonic approximation to a v is natural at higher 
densities, and could be quite good at low T, 
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a v ((f),T,p) = '-j-^^ln^haji))^) 

i=l 

, 3N 

i -M^o) + M^M>))(£)], (4) 
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where Ui is the i'th normal mode frequency, ujq is the 
frequency unit, and the averages are over representative 
IS with energy <p. Sastry 0] found that the frequency 
sum has a linear </>-dependence, 



f 



3JV 

- X J2 in(wi/«o))(0) = f°(p) + f l {p)ct>, (5) 

i=l 



for some liquid states of the Lennard-Jones (LJ) mixture. 

La Nave et al. Q reached the same conclusion for OTP, 
and obtained and tested the Gaussian - linear harmonic 
(denoted GLH) equation of state. With these approxi- 
mations, 



<j>(T,p) = (Mp)-f(p)5 2 (p)) 



S 2 (P) 
k R T' 



(6) 
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The pressure has the form 



P(T,p) = TP T (p) + P const (p) + T- l P 1/T (p), (7) 



where 



V^(Sc,oc-M/° + /Voo)) 
Pconst = p 2 {d(j>oo/dp) 



. d , S 2 . 



I'r = -/>-tH-v.. - I) (8) 

(9) 
(10) 

The 75 energy reaches a high-T plateau, = 0o — / 1( 5 2 , 
and, correspondingly, s CiOC = k B o-(4> QO , p). 

The "potential energy landscape equation of state" , 
evaluated via computer simulation, was shown to accu- 
rately reproduce the true pressure for a range of liquid- 
state p and T in OTP. Here we attempt to reach some 
more general conclusions. We identify the origin of the 
van der Waals equation in the IS formalism, and ob- 
tain expressions for density-dependent "a" and "6" co- 
efficients. There is an additional contribution, P\/t in 
the GLH approximation, which we call the "landscape 
pressure". The resulting extended van der Waals equa- 
tion has the possibility of yielding a water-like density 
anomaly (already pointed out by Sciortino et al. |ll|L 
flat isotherms in the coexistence region, as opposed to 
vdW loops, a positive derivative {dU/dp)T at high den- 
sity, and a critical anomaly in Cy . We make use of the 
GLH approximation but also derive results independent 
of any assumptions about the vibrational free energy. 
Relevant PEL quantities are determined by computer 
simulations on the single-component LJ fluid for a wide 
range of densities including the coexistence region. Fits 
to the data lead to to an analytic equation of state. 



II. THE EXTENDED VAN DER WAALS 
EQUATION OF STATE 



In the vdW equation the pressure, 

k B Tp 



P(T, P ) = 



1 - bp 



ap 



(11) 



is a sum of a T-independent term and a linear term. Here 
a and b are the well known coefficients expressing the 
reduction of P from the ideal gas value by the attractive 
forces and the increase due to repulsions, respectively. 

1. The Gaussian linear harmonic (GLH) approximation 

Identification of the attractive and repulsive-ideal gas 
contributions with P CO nst and Pt is obvious, and holds 
up upon further consideration. The standard expres- 
sion is a = —(d(U/N)/dp)T, where U is the total en- 
ergy. In the GLH approximation, U — Ui S + ^NkgT, 
and (d(JJ / N) / d p)t — (d<j)/dp)T- In general the deriva- 
tive is T-dependent but the plateau value, {d4>oo/dp)T, 



is a function of density only. Thus (Eq|5J) we propose a 
density-dependent extended vdW a coefficient 



Pconst = ~a(p)p 2 

a (p) = -(d(poo/dp) 



(12) 
(13) 



and the less familiar IS energy is related to a textbook 
parameter. 

The vdW repulsive/ideal gas pressure is conventionally 
derived from the entropy. In the GLH approximation it 
is easy to see that Eq[S]just involves the derivative of the 
total entropy/particle, s = s c + s v , and we propose 



Pi 



20/ 



pki 



i-pb(py 



(14) 



with a density-dependent b coefficient. Since the entropy 
is evaluated on the plateau, Pt is T-independent. The 
contribution to the pressure is TPt, with the explicit 
multiplicative factor coming from the way entropy enters 
the free energy, —TS. It might seem that harmonic oscil- 
lations have nothing to do with the molecular congestion 
that increases the pressure at high density, but the den- 
sity dependence of s v is entirely due to the change in the 
shape of the basins, expressed through the frequencies, 
and that is relevant to packing. 

Eq EH does not look terribly transparent. For clarifi- 
cation consider that, if 



Soo = k B ln{-f- — ), 

P 



(15) 



where 7 is a p-independent constant, and if for now we 
let b be p-independent as well, the vdW Pt is obtained. 
Eci ll5l is plausible, with entropy/particle properly varying 
as ln(l/ p) at low density and vanishing as "close packing" 
is approached, p — > 1/6. Specifically, s M vanishes at 
p = (b + I/7) However real molecules with soft cores 
retain positive entropy as they are compressed to very 
high density, so any vanishing of S must be understood 
as an extrapolation from a particular density range. We 
anticipate that b{p) defined in EqE3 w nl decrease at high 
p, keeping the entropy positive and the pressure finite. 
This expectation will be explicitly realized in Sec II VI 

There is considerable interest [3 IT3I [l4| in the varia- 
tion of the parameter a, which determines the total num- 
ber of IS, with density and from substance to substance. 
If the p-dependence of a (via <t(0oo)) dominates that of 



1 / 1 - Pb, 
a w m(7 ). 

P 



(16) 



Thus, the free energy evaluated on the high-T plateau 
of the PEL yields an equation of state with more com- 
plicated density dependence than the vd W equation, but 
the same T-dependence. The extended version also takes 
into account the deviation of <j)(T) from its plateau value, 
measured by d(j)(T,p) = <poc{p) ~ (f>{T,p), In the GLH 
approximation, where d<fi(T) — S 2 /ksT, this gives rise 
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to Put oc (dS 2 /dp). With no density dependent vibra- 
tional energy d(f>(T) is the entire relevant energy drop, 
and it also controls (S Ct00 — S C (T)) via a{4>oo — d(f>(T)); 
everything is proportional to 5 2 . Adding the landscape 
pressure yields the extended vdW equation in the GLH 
approximation, 



P(T, P ) = 



k B Tp 
1 - b(p)p 



d 



-a{p)p 2 -p 2 -{ 



dp y 2k B T 



) (17) 



2. Beyond the GLH approximation 

Simulation shows that the GLH expression for 4>(T) is 
only a crude representation; what is the source of the er- 
ror? The Gaussian approximation has no upper or lower 
bound on possible IS energies, but is a good representa- 
tion for the contributing fluid-like states over a substan- 
tial range of T and p, and has 0, Q, E| some theoretical 
justification. 

The harmonic approximation is valid for a particular 
normal coordinate when the system remains close to the 
IS. This suggests high density, where in fact the harmonic 
4>(T) is better, and deterioration as p is decreased to- 
wards the coexistence region. The imaginary frequency 
instantaneous normal modes |l5l |. whose number corre- 
lates [H E [H El very well with the self-diffusion 
coefficient D and increases with decreasing density, are 
completely anharmonic vibrational coordinates. Thus we 
aim to retain a Gaussian density of states, but to aban- 
don the linear harmonic approximation for a v . Deriving 
specific expressions for a v is a current focus of our re- 
search but we will now just try to draw some simple, 
general conclusions. 

The true <j>(T) does reach a high-T plateau. Accord- 
ingly, we rewrite the free energy, 



A(T,p)/N = [c/ioo + Oc^oo, T,/o) - k B Ta(^oo,p)] 



-[#(!» - k B Td0 2 (T,p)/25 2 } - [a. 



,T,p) 



a v (cj)(T, p), T, p) - k B T(4> - oo )d^i(T, p)/S 2 



(18) 



suppressing the p-dependence of the constants 0oo, 4>q 
and 5 2 . EqElhas three distinct groups of terms, en- 
closed by square brackets. The first group is the plateau 
free energy, which yields P CO nst and Pt in the GLH ap- 
proximation, the second group describes the effect of the 
drop of 4>(T) from the plateau and becomes the landscape 
pressure Put in the GLH, and the third group vanishes 
in the GLH. 

As a first step consider the form 



a v (<t>,T,p) = (e° v (p) + el(p)<P) 

-T(s° v (p) + sl(p)<l>) + f(Tonly); 



then. 



<P(T) = (0o 



,1X2 



-) 



(l+el)6 2 



(19) 



(20) 



Eq^]is hardly the most general possibility but it does 
include linear ^-dependence of both the vibrational en- 
ergy and entropy, while the linear harmonic approxima- 
tion has nontrivial entropy only (f 1 — —s\Jk B , Eq^Jl. 



Of course, a function of T only does not influence P or 
4>(T). With the Gaussian, Eall9lconstitutes the Gaussian 
linear (GL) approximation. 

Returning to Eg 1181 the first group of terms now con- 
tains the total (configurational plus vibrational) energy 
and entropy on the plateau. The vibrational entropy 
terms in the third group still cancel but the deviation 
de v (T) of the vibrational energy from the plateau value 
survives, and combines with d<fi(T, p) in the second group 
to form the deviation of the total energy from the plateau 
value. The new contribution de v (T) will vanish as the 
harmonic approximation becomes accurate, i.e, high den- 
sity. Evaluating the second group with Eq [5D| the In- 
dependence of P is unchanged from the GLH approxi- 
mation and 



p „' 

1 const — P 



dp ' 



Pi 



2 d 



P o (^c.oo Sv.oo} 

dp 



Put = -P'jri 



2 d Al + el) 2 5 2 



)• 



The corresponding GL extended vd W equation is 



P(T, P ) 



k B Tp 



1 - b(p)p 
Ea [HI still holds for b, but 



d Al + el) 2 5 2 



(21) 
(22) 
(23) 

■)• (24) 



a{p) 



dp 



dp 



<V 2fc B T 



((l + ei^oo+e"). (25) 



A more general treatment of anharmonicity is neces- 
sary to discuss densities, including the liquid-gas coexis- 
tence range, where <f>(T) is not described by Eq[20| In the 
GLH and GL approximations there is no characteristic 
temperature at which the plateau is reached; d<p(T) van- 
ishes smoothly as 1/T. Simulation, however, shows lid 
an identifiable "top of the landscape" temperature |19| , 
TtoL-, at typical liquid densities. Considering the exact 
Eq^Jfor the free energy, if we assume that true plateau 
quantities are available, they will yield P con st and Pt ac- 
cording to Eqs |^ and [23 It is not terribly difficult to 
obtain the true <j>(T) from computer simulation, so let us 
also consider that it is known; then the second group can 
be evaluated. The third group, with its cancellations, is 
less certain. In the GL approximation it converted d(f>(T) 
to d<j)(T) + de v (T). If we assume that happens more gen- 
erally we obtain a third Gaussian-anharmonic extended 
vdW equation of state 



P(T,p) = -^--a(p)p 2 - 
1 - b{p)p 

k B Tdtf(T,p) 



p z — (d4>{T,p)+de v {T 1 p)~ 



2S 2 



). (26) 



The coefficients a and b are given by Eas 1251 and 1141 re- 
spectively, with the understanding that the true thermo- 
dynamic functions are to be employed. Their properties, 
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and the behavior of </>(T), implicitly incorporate the ef- 
fect of an anharmonic a v upon the pressure. 

Again, the argument is that, since the GLH approxi- 
mation for dcj)(T) can be quite poor, the most important 
step is to correct the contributions to P which directly 
reflect d<j)(T). On the other hand, contributions evalu- 
ated on the plateau are left in the extended vdW form. 
They also contain anharmonic effects implicitly through 
any non-GLH behavior of (j)oo(p) and s((/>oo, p) when they 
are evaluated, e.g, via simulation. 

III. SOME GENERAL PROPERTIES OF THE 
EXTENDED VAN DER WAALS EQUATION 

The extended vd W equation has density-dependent a 
and b coefficients, but its most striking feature is the 
landscape pressure, in which we now include the de v (T) 
contribution, with a new T-dependence. This term is 
oc 1/T in the GLH and GL approximations, has a more 
general form in Eq|22| and physically expresses the ef- 
fect on the pressure of the descent of the system from 
the high-T plateau in the IS energy. Thus we have a 
PEL-b&sed criterion: the extended vd W equation should 
provide a good representation of thermodynamic states 
corresponding to the plateau. In the Angell classifica- 
tion of liquids on a strong- fragile scale [lj| 1211 l2lL |22| , 
the more complex behavior of fragile liquids is associated 
with a more pronounced drop d<j){T) from the plateau. 
Thus we further suggest that the ordinary vdW equation 
should be most useful for strong liquids. 

1. Isotherms 

Perhaps the most-discussed features of the vdW equa- 
tion are the loops in the isotherms below the critical 
temperature T c . They are a property of the metastable 
or unstable homogeneous fluid in the coexistence region. 
There is no provision in the vd W equation for the phase 
separation which goes hand in hand with the true, flat 
isotherms. 

There are two possibilities upon crossing the liquid- 
gas coexistence curve, Ti g (p), from above in the (p,T) 
plane: the system may phase separate in thermal equi- 
librium, or it may remain a metastable, homogeneous 
fluid. Similarly, above the triple point density, cooling 
below the melting temperature may lead to liquid-solid 
phase separation or a metastable supercooled liquid. At 
a given density, there will exist IS representative of all 
the possible thermodynamic states. Under the coexis- 
tence curve, there will be both homogeneous fluid and 
phase-separated liquid-gas and gas-solid IS. Above the 
triple point density the IS types will be homogeneous liq- 
uid, liquid-solid phase separated, and crystal with vary- 
ing amounts of disorder. Due to surface effects, phase- 
separated IS may be difficult to observe, and/or modified 
in character, in finite- N simulation. A FiSL-based calcu- 
lation in a metastable state is achieved by including only 
the IS to which the system is restricted, while equilibrium 
results from keeping all the IS. 
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FIG. 1: Ordinary and extended vd W pressure vs density for 
(T c -T)=0.015, a=4, 6=1 (p c =l/3); all quantities in LJ units. 
Outside the coexistence region the two pressures are identical, 
inside the flatter curve is the extended vdW equation with the 
illustrative form of d<j>(T, p) from the text. 



Phase separation causes a large drop in <fi(T), or a 
strong increase in d(/)(T). On the other hand d(j>(T) 
varies more gently when the system remains homoge- 
neous. Thus the landscape pressure in Eq 1261 behaves 
quite differently for the two cases. We suggest that it is 
small for homogeneous states, making the extended vdW 
equation a good approximation, but is large and leads to 
flat isotherms if phase separated states are included in 
the evaluation of d<f>(T). 

To demonstrate this idea, let phase separation at a 
given p, and the corresponding strong growth in dcf)(T), 
begin at T/ a (p). For T < T c , the distance to the coexis- 
tence curve, (Tig(p) —T), is an increasing function of p 
for p < p c and a decreasing function for p > p c ; the den- 
sity derivative changes from positive to negative at p c . If 
d<p(T) is a monotonic function of (Ti g (p) — T), its deriva- 
tive (dd<p(T)/dp) will have the same behavior. Then the 
landscape pressure will be negative below p c and positive 
above, i.e. it will form a van der Waals anti-loop which 
can cancel the vdW loop and produce a flat isotherm. 

As an illustration only, consider the simple anzatz 
d<j)(T,p) = c(Ti g {p) - T) 2 , with Ti g estimated from the 
ordinary vd W equation. We use p-independent a=4 and 
6=1 for an approximate description of the LJ fluid and, 
since the growth of d(f>(T) is the physical effect of inter- 
est, ignore de v . All quantities are expressed in natural LJ 
units. The landscape pressure is nonzero within the vdW 
coexistence curve only, so we assume it does not change 
the critical point from T c = 1.185, p c = 0.333. Fig Q] 
shows vdW and extended vd W isotherms at T = 1.170 
for the choices c=1.91, S 2 =l. We do not claim to have the 
correct d<p(T) but the point is that any model in which 
its growth begins at the coexistence curve will yield an 
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anti-loop and potentially a flat isotherm. 

2. Some thermodynamic derivatives 

There is considerable current interest in the phe- 
nomenon of a negative (dV/ dT) p, the "density anomaly" 
well known in water, which may or may not be associated 
|23| with the existence of multiple critical points. Some 
textbook manipulations show that equivalent conditions 
are {dS/dV) T < or (dP/dT) v < 0. 

It is immediately apparent that the vdW equation can- 
not have a density anomaly, since 



(dP/dT) v 



k B p 
I -bp 



> (van der Waals), (27) 



but things are quite different [Tlj when the landscape 
pressure is added. In the GL approxination, 



(dP/dT), 



k B p 



e\?8 2 



1 - b(p)p 



1 dp K 2k B T 2 



(28) 



and thus, if the density derivative is negative and large 
enough, thwre will be an anomaly. A more complicated 
condition can be expressed with the true d<p(T) via Ea l26l 

There exists a clear physical interpretation of why the 
extended vd Inequation can have a density anomaly. One 
expects that entropy should decrease with increasing den- 
sity (positive (dS '/ '8V)t), as the system becomes more 
congested, but this is only true on the high-T plateau 
where S is a maximum. In the GL approximation, the 
deviation of <j>(T) from the plateau at constant T is pro- 
portional to S 2 . If 5 2 decreases with increasing density, 
the system gets closer to the maximum plateau entropy, 
which may compensate for the decrease in Soo itself, lead- 
ing to (dS/dp) T > and (dS/dV) T < 0. 

For the vdW equation the energy decreases with 
increasing density; the relation (d(U /N)/dp)x = 
T(dP/dT) v -P yields {d(U/N)/dp) T = -a, correspond- 
ing to the negative energy of attraction. This is because 
repulsions enter via the entropy only, i.e, there is no true 
positive repulsive energy contribution to the pressure, no 
matter how high the density. In the extended vd W equa- 
tion, (dejjoo/dp) can fSec IIV|I become large and positive, 
corresponding to a large negative a(p), at high density; 
the "vdW attractive" term then transforms into a repul- 
sive pressure. Furthermore, taking the landscape pres- 
sure into account, in the GL approximation, 



(d(U/N)/dp) T = -a(p) 



Aril 



-), (29) 



which can possibly change sign for the same reason given 
in the last paragraph for the density anomaly. 

The critical behavior of the heat capacity is described 
[24l | by the exponents a and a', 

(dU/dT) Vc oc (T-T c )- Q , T>T C 
cx(T c -T)- Q ', T <T C 

In the vd W approximation |24j, a=a'= 0. Since we are 
obtaining P from the free energy the simplest consistent 




FIG. 2: Simulated Ttol(p) data (jagged) and coexistence 
curve from Ref 251, LJ units. 



route to U is also through A; elementary statistical me- 
chanics yields the exact relation 



(U(T)/N) = cf>(T) 



2 d(a v (T,cf>)/T) 
[ dT 



(30) 



where the notation is intended to convey that the deriva- 
tive acts on the explicit T-dependence of a v only, not on 
that which enters implicitly through 4>(T). For an equa- 
tion of state based upon a particular approximation to 
a v , use of that approximation in Eg I3UI yields the cor- 
responding potential energy. Since the IF equation is 
obtained by ignoring the deviations of <f> and e v from 
their plateau values in a v , strong T-dependence of U(T) 
is discarded. That is why a=a'— 0; the critical anomaly 
is in the landscape term. 



IV. APPLICATION TO THE LENNARD-JONES 
LIQUID 

Here we present some preliminary results on LJ. A 
careful study of phase equilibria [2j| determined that the 
critical point and triple point are p c =0.31, T c =1.31 and 
p t =0.84, Tj=0.75, LJ units. An estimate of Tt l(p), de- 
termined from 4>(T) , N=256, is shown in Fig |2 along 
with the coexistence curve from Ref |25t | . A minimum is 
found near the triple point density, and with decreas- 
ing density Tt l(p) runs between the liquid-gas spin- 
odal and coexistence curve to the critical point. Earlier 
we discussed the liquid-gas coexistence temperature, Ti g , 
as marking the onset of a rapid rise in d<f>(T); Ti g is a 
" ToV temperature for equilibrium states. Practically, 
the metastable homogeneous phase can survive until the 
spinodal is reached, and our result is expected in finite- 
size simulation. What may be interesting is that the 
liquid-gas spinodal and the ToL temperature discussed 
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FIG. 3: Plateau energy (j>oo (lower), shifted up by 5 energy 
units, and squared Gaussian width S 2 vs reduced density, all 
in LJ units 



in supercooled liquids, found at densities above the min- 
imum, are thus connected. Leyvraz and Klein |27| have 
suggested that properties of supercooled liquids may be 
influenced by a spinodal. 

We have obtained S 2 (p) and </>oo(p) (Fig El f° r 
2.0>p>0.20, N=500, from the distribution of IS visited 
at high T, where 5 2 becomes T-independent, and from 
</>(T), respectively. Starting at low density both quan- 
tities, like TtoL, decrease to minmium values near the 
triple point, and then begin to rise. Similar behavior has 
been observed 0, Hi| for the IS pressure. Correspond- 
ingly the vdW attractive term becomes a positive, repul- 
sive pressure. In the GLH approximation the behavior of 
S 2 (p) would cause the landscape pressure to vanish near 
the triple point, identifying a region well described by 
the extended vdW equation. 

The higher the density, the better we expect the har- 
monic approximation to perform. Thus we are going to 
use our data to evaluate the GLH equation of state over 
the entire available density range, but but we do not make 
any claims of validity at low to intermediate density. Fit- 
ting S 2 (p) and (/?) to a sum of two exponentials gives 
good fits, particularly at high density, 2.0>p>1.2. Taking 
density derivatives yields P CO nst and the landscape 

pressure P\/t- 

We do not at the moment have the information nec- 
essary to calculate b(p) from first principles. This coef- 
ficient involves an extrapolation, i.e, in some small den- 
sity range it appears that the pressure would diverge at 
a particular, higher "close packed" density, but since the 
cores are soft that density is never reached. Consequently 
b may appear to be constant at low density, but its p- 
dependence becomes essential at high density. 

It is difficult to equilibrate the system at p>1.2 without 
going to quite high T. To get some idea of the behav- 
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FIG. 4: High density b coefficient from simulation and fit to 
EaEiTl &o=1.31. Note 6(/9 c =0.31)=1.08. 

ior of b we have calculated the T=25 isotherm, and b(p) 
from the GLH Eq El The results are shown in Fig 0] 
Because of the soft cores, the quantity (1 — b(p)p) can 
become small, but never zero. Thus we suggest a plausi- 
ble behavior is exponential decay with p, giving 

b(p) = (1 - e- b °»)/p, (31) 

and the smooth curve in Figgis Eq|^with &o=1.31. 
Combining everything, 

P(T, p) = Tpe 1 - 31 ? + p 2 (9.15e 195p - l3.7e 1A8p ) 

_£_( .0168e 4 - 96 " - 35.4e- 3 ' 09 ), (32) 

anticipating that the analog of the vd W attractive term 
will be positive, and the landscape term negative. 

Matching the [2j| true p c and T c of LJ to the vd W equa- 
tion gives b—1.07 and a=4.73 At the critical density we 
find a(0.31)=4.92 (negative of first parenthesis in Ea l32|) 
and (Ea l31|) b(p)=1.08. Again, we expect the missing an- 
harmonicity to be important at the critical density and 
regard these results with some skepticism, but the agree- 
ment is remarkable. The density dependence of a, along 
with b, is shown for p < 0.90 in Fig |3J note the sign 
change near the triple point density. 

V. SUMMARY 

The equation of state is obtained from the density 
derivative of the Hclmholtz free energy, as always. The 
PEL approach simply provides a less traditional way to 
view the problem and generate approximations. Since a 
Gaussian approximation for the IS energy distribution is 
reasonable, the focus is on the vibrational free energy a v . 
The GLH approximation was worked out 5] by La Nave 
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FIG. 5: a (upper) and b coefficients vs density 

et al., and the starting point of this paper is simply the 
observation that their result resembles the van der Waals 
equation with density-dependent a and b coefficients, and 
an extra "landscape" term. 

In addition, we have tried to include the anharmonic- 
ity of a v with as few assumptions as is possible; in the 
extended vdW Eg 1261 anharmonicity is implicit in the In- 
dependent inherent structure energy, and in the plateau 
entropy. The smaller the drop of the IS energy from its 
high-T plateau, the smaller the landscape pressure, sug- 



gesting that the vdW equation is best suited to strong 
liquids. The extended equation can reproduce, and pro- 
vide PEL interpretations of, thermodynamic phenomena 
absent from the usual vdW equation. 

Computer simulation in the LJ fluid yields a(p) and 
b(p). The a coefficient becomes negative at p > 1.0, as 
the "vdW attractive" pressure becomes a repulsive pres- 
sure. The b coefficient is represented by a simple ex- 
pression stemming from the idea that the pressure may 
become exponentially large, but not infinite. Since the 
harmonic approximation may be accurate at high den- 
sity, the resulting analytic expression is suggested as a 
high density equation of state, which we will explore in 
future work. Nevertheless, when evaluated at the criti- 
cal density, a(p c ) and b(p c ) are remarkably close to the 
accepted values. 

The quantities 0oo (p) and S 2 (p) have minima near the 
triple point. Thus in the GLH approximation both P con st 
(a) and P±/t (landscape pressure) vanish. The "top of 
the landscape" temperature also has a minimum near the 
triple point. At lower density Tt l(p) runs under the co- 
existence curve to the critical point, at hi gher density it 
rises and becomes the quantity discussed .19] for super- 
cooled liquids; this connection may support the conjec- 
ture |2jj of Klein and Leyvraz. 
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